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Abstract 



The effective dc-conductivity problem of isotropic, two- 
dimensional (2D), three-component, symmetric, regular com- 
posites is considered. A simple cubic equation with one free 
parameter for a e {ffi, oi, 03) is suggested whose solutions au- 
tomatically have all the exactly known properties of that 
function. Numerical calculations on four different symmet- 
ric, isotropic, 2D, three-component, regular structures show 
a non- universal behavior of a a {o\, (72, C3) with an essential 
dependence on micro-structural details, in contrast with the 
analogous two-component problem. The applicability of the 
cubic equation to these structures is discussed. An exten- 
sion of that equation to the description of other types of 2D 
three-component structures is suggested, including the case 
of random structures. 

Pacs: 72.15.Eb, 72.80.Tm, 61.50.Ah 



I. INTRODUCTION 

The classical duality transformation of two-dimensio- 
nal (2D) heterogeneous composites, discovered by Keller 
jjj and independently by Dykhne j| , has been applied in 
a restricted set of physical contexts. The dual symmetry 
is based upon the observation that any 2D divergence- free 
field, when rotated locally at each point by 90°, becomes 
curl-free, and vice versa. This leads to the result that 
static effective physical properties of 2D infinite hetero- 
geneous composites, like electric conductivity a e , ther- 
mal conductivity, dielectric permittivity, as well as some 
other static properties, satisfy some exact relationships 
which follow from the similarity between dual problems. 
For example, in the case of a two-component composite, 
a universal square-root law behavior of the bulk effective 
transport characteristics was demonstrated 0, ||. It is 
also known[] that this transformation leads to exact re- 
sults in 2D infinite composites made of an arbitrary num- 
ber of components. Dykhne 0] gave sufficient conditions 



1 It seems to be strange but we have not found throughout 
the papers concerned with this subject any published proof 
of this statement. Such proof is so useful that we give it in 
Appendix. 



that, when satisfied by the component conductivities <7j, 
lead to exact results for the isotropic bulk effective con- 
ductivity o- e (ai). We know of only one attempt Q to 
consider a 3-component 2D composite with a doubly pe- 
riodic arrangement of two kinds of circular inclusions em- 
bedded into the matrix. The effective permittivity was 
obtained using the dipole approximation for the inclusion 
polarizations. But we do not know any rigorous results 
obtained for any 2D composite microstructure, made of 
3 (or more) components, with arbitrary component con- 
ductivities. Apparently, this is not an accident but re- 
flects the disappearance of commutativity in symmetric 
groups when we upgrade from the ^2 permutation group 
to the S n permutation group. 

In this paper we consider the effective isotropic con- 
ductivity problem for 3-component 2D infinite compos- 
ites with translational order, with a microstructure that 
is symmetric in the 3 components. We formulate an 
approach based on the conjecture of the algebraicity of 
a e (ai, (72, 03) and its general properties. We have found 
that the algebraic equation of minimal order where all 
these properties can be satisfied is a cubic equation, 
which contains 1 free parameter, with coefficients made 
of the independent invariants of S3. This equation is in 
agreement with Dykhne' result (see Eq. ([?])). Its pre- 
dictions are compared with numerical solutions for a e 
in some regular 3-component microstructures. It is also 
extended so as to apply to other types of 3-component 
composites, including random microstructures. 



II. TWO-DIMENSIONAL THREE-COLOR 
COMPOSITES 

The effective dc-conductivity problem for n-component 
symmetric, 2D, infinite composites with translational or- 
der, can be reformulated with the help of n- color plane 
groups. 

Color groups are generalizations of the classical crys- 
tallographic groups. Different colors may correspond, for 
example, to different chemical species or, more generally, 
to different values of a physical property which is defined 
as a tensor of /c-rank : scalar - density p , vector - mag- 
netic spin m , tensor of 2-nd rank - conductivity Uij, 
tensor of 3-rank - piezoelectric modulus dij n , etc. 
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Every n-color plane group has its origin in one of the 
Ni =17 regular (color-blind) plane groups. The number 
N n of n-color plane groups is a non-monotonic function 
of n: N 2 = 46 §; N 3 = 23; N 4 = 96; N 5 = 14; N 6 = 90 
H , • The plane groups n < 60 are tabulated in S . 

Only 10 of the 23 3- color plane groups have a 3-fold 
rotation axis: 

5 lattice equivalent groups - P3(L)|P1(L), P6(L)|P2(L), 
P31m(L)|Cm(L), P3ml(L)|Cm(L), P6mm(L)|Cmm(L) and 
5 class equivalent groups - P3(L)|P3(L ), P6(L)|P6(L'), 
P31m(L)|P3ml(L'),P3ml(L)|P31m(L'), 
P6mm(L)|P6mm(L ) where L is a possible sub-lattice^ of 
L invariant under rotations of order 3 §. All these groups 
are compatible only with hexagonal Bravais lattices. The 
3-fold rotational symmetry makes the effective conduc- 
tivity in structures governed by those groups isotropic. 
This follows from the Hermann theorem [0] about fc-rank 
tensors in media with an inner symmetry which includes 
a rotation axis of highest order r, r > k. Despite their 
different geometries, all these structures have one impor- 
tant property in common: They are all invariant under 
the full permutation group S3 which exchanges the colors, 
therefore they are related to the permutational crystal- 



lographic color groups [10 



III. EFFECTIVE CONDUCTIVITY a R (a x , <7i, ffs) 
AND ITS ALGEBRAIC PROPERTIES 

A direct way to solve the dc-conductivity problem for 
an n-component composite begins with the local field 
equations 

V x E(r) = , V ■ J(r) = , J(r) = cr(r) E(r) , (1) 

along with appropriate boundary conditions for the elec- 
trical potential. The local conductivity cr(r) is a discon- 
tinuous function er(r) = a i} if r € Aj , i = 1,2,3 where 
Ai is a homogeneous part of the composite with constant 
conductivity er^. The isotropic effective conductivity o~ e 
can be defined via Ohm's law for the current J e and the 
field E e averaged over the system 

J e = a e E e , J e = i J 3(r)dS , E e = i J E(r)dS . 

(2) 



2 We follow the notations of Ref. jj| for three-color plane 
groups Gcoi = G(L)|G (L ) which means that G is the geomet- 
rical, or Fedorov, plane group and the subgroup G C G of 
index 3 contains operations that keep the first color fixed. 
We do not specify here the relationship between the lattice L 
and its sub-lattice L . 



Except for a medium with trivial ID inhomogeneities (a 
layered medium) an exact solution of this problem does 
not exist for any regular or random structure. At the 
same time the function o e (pi) must have the following 
general properties which we are going to exploit: 
1. Homogeneity of 1 — st order 



a e [k oi, k o- 2 , k er 3 ) = k a e [a x , er 2 , 03) 



(3) 



This follows from the linearity of the static Maxwell equa- 
tions (0) and from the definitions of the average current 
and field (|). 

2. Permutation invariance 



o- e (Pl {0-1, 0-2, 0-3}) = o- e (a l7 er 2 , a 3 ) 



(4) 



where Vi, I = 1,...,6 is a permutation operator of the 
indices {1,2,3} (or of 3 colors). The six operators Vi 
form the non-commutative group S3. The existence of 
permutation invariance presumes that the 3 components 
are distributed with equal volume fractions pi = p 2 = 

= 1/3. 
3. Duality 



c e (eri, c 2 , 0-3) X er e ( 



1 1 



1 



= 1 



U\ cr 2 CT 3 



SeeEq. ( [All ). 
4. Compatibility 



a e (a, a, a) = a 



(5) 



(0) 



The last formula does not follow from the previous ones 
but reflects a natural requirement. 

Dykhne Q proved a theorem for a symmetric compos- 
ite of 3 components, namely 



cr e (ai, er 2 , 03) = \J0\02 , if 



(7) 



The last formula represents the only known rigorous re- 
sult for 3- color 2D isotropic composites. It can easily be 
shown to follow from (^) and (J5J) . Let us mention one 
more conclusion which follows from (^) 



a e (a, 0, 0) = , 



(8) 



which reflects the percolation property of such a compos- 
ite. In fact, it is easy to construct 3- color 2D isotropic 
composites with the more specialized property 



a e (a 1 , cr 2 , 0) = 



(9) 



Fig.[j]d illustrates such a structure with traps, i.e. a net- 
work of closed, simply connected loops, which enclose 
plaquettes of the two other colors. If any color denotes 
an insulator, then the effective conductivity of the com- 
posite is clearly zero. The 2- color 2D square checkerboard 
has a similar property. It differs from the case of Fig.gd 
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by the value of the percolation threshold. The previ- 
ous example (|9|) shows a non-universal behavior of a e (<Ji) 
even for symmetric isotropic microstructures, with essen- 
tial dependence on micro-structural details, in contrast 
to the 2-color 2D symmetric isotropic composites, where 
o e = ^]<J\Vi always || . 

We will look for o e (o%, o%, 03) satisfying the require- 
ments (|3|- [|) among algebraic functions. This choice is in- 
spired by the fact that the symmetric 2-color 2D isotropic 
composite generates a quadratic equation for o e (o\, 02). 
Another motivation is the situation which exists vis-a-vis 
some discrete 2D models in statistical mechanics, where 
a duality transformation exists that is similar to the one 
invoked here. Such a transformation exists for the Ising 
and Potts models |ll]] , where the critical point equations 
are algebraic — quadratic equation for the Ising model 
on a self-dual square lattice and cubic equations for the 
Potts model on the mutually dual triangle and honey- 
comb lattices. 

According to the "fundamental theorem of symmet- 
ric functions" the symmetric group S3 has 3 alge- 
braically independent homogeneous invariants (basic in- 
variants) 



il = Ol + (7 2 + OS , 1 
I3 = 0\ (72 (J 3 , 

which satisfy the obvious restrictions 



(71(72 + (72(73 + (73(7i 



(10) 



1 h 
-h > — 

3 " h 



> 3 



1 



> 3 



and can be used as independent variables instead of 
01 j °2 j 0-3 . The difference I2 — h \fh can have either 
sign. The function o e (I\, I2, h) now satisfies (^) auto- 
matically. 

It can be shown that the algebraic equation of mini- 
mal order for the function o e (Ii, I2, 13), which satisfies 
all basic requirements, is a cubic equation of the form 



A Ii o e — A I 2 o e 







(11) 



where A is a free parameter responsible for the non- 
universality, and <7 e is a value bounded from both above 
and below (ljj 



>3^ 
h 



(12) 



It is easy to show that the equation (|llj) satisfies all basic 
requirements (||^) and automatically satisfies Dykhne's 
theorem (7) independently of A. Indeed, the require- 
ments (||, I, ||) one can check immediately. The duality 
property (5) follows when we notice that (11) is equiva- 
lent to the cubic equation for — ( — , — , — ): 



1 4/1 1 1 ^ 

— + A(— + — + —) 

0% (7l (7 2 CF3 



A(- 



1 



1 



= 



(7l(72 (72(7 3 (7 3 (7l' 



(7l(72(7 3 



The property (]?]) can be proven by straightforward sub- 
stitution (73 = ^/(7i(72 into ( |TT| ) . 

As we will see later, A reflects not only the plane group 
of a color tessellation but also the shape of the elementary 
cell. To illuminate the meaning of A let us put into ( p"l| ) 
A = |. After simple algebra one obtains the following 
equation 



0~e - 0~1 Q~e — Q~2 °~ e ~ & 3 
(7 e +(7i (7 e +(7 2 °~e+°~3 







(13) 



which coincides with the Bruggeman effective medium 
approximation for a symmetric, 3-component, 2D com- 
posite fuj. 

From restrictions ( |l2] ) and after straightforward ma- 
nipulations, one can show that A is bounded from below 



oo > A > 



(14) 



It is noteworthy that the lower bound ensures that (|ll| ) 
has only one positive root, thus avoiding the possibility 
of multiple physical solutions. More accurate estimation 
of the lower bound for A, using the Hashin-Shtrikman 
|]l5| exact bounds cr^ s for isotropic conductivity o~ e of 3- 
component 2D composite, does not change (|l4]). Indeed, 
for (7i > (72 > (73 one has 



h > o-Jj S > (7 e > cr HS > 3 



h 
h 



(15) 



where 



'HS 



o-i ■ (1 

t hs = o-3 • (1 



4(cr 2 o- 3 - of) 



(72(73 + 3(71 ((72 + O3) + 5of 

4(0102 - ol) 



<7l(72 + 3(73 ((71 + (7 2 ) + 5cr| 



)• 



Due to the ambiguity of the general case, we consider the 
special case o\ > 0% = 03. Here we have 



O hs = (71 



5 + r 
1 + 5r 



O hs = (7 3 



l + 2r 



(71 

— > 1 

O3 



We can then derive the following lower bounds from (|llj) , 



A > 



l + 5r 1 



A> - 



l + 2r (2 + r) 2 : 
r 2 + r 1 



18r 



1 + 5r 1 + r 1 + 16r + r 2 



Since the last inequalities must fit all cases including r 
00, we return to (111). 



IV. EFFECTIVE CONDUCTIVITY OF REGULAR 
STRUCTURES: NUMERICAL RESULTS 

In the present section we will study numerically four 
different infinite 2D 3- color class equivalent regular struc- 
tures of P6mm(L)|P6mm(L') and P6(L)|P6(L') types. In 
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order to avoid cumbersome notations for these struc- 
tures, we will use the following symbols: He (honeycomb) 
(Fig.0a) ; Fl (flower) (Figj^b) ; Co (cogrose) (Fig.0c) ; Rh 
(rhombus) (Fig{j]d). 

Due to the homogeneity (^) of Eq. ( pi] ) one can rescale 
<7 e (ai, (72, 03) defining a new function 4(4>4) where 
(j- = <Ji/oi, i = e, 2,3. A typical shape of the surface 
defined by 4(4; 4) is shown in Fig. 0. 




(a) (b) 




(c) (d) 



FIG. 1. 3-color 2D composites: a) He 

-G col = P6mm(L)|P6mm(L'), b) Fl -G col = P6(L)|P6(L), c) 
Co -G co i = P6(L)|P6(L), d) Rh -G co , = P6mm(L)|P6mm(L'). 




FIG. 2. Typical shape of the surface a e (<j' 2 , a'-s) for 3-color 
2D isotropic composites. The loop where Dyktme's theorem 
is satisfied is shown by the dashed lines. 



This surface always contains the following points 
4(0,0) = 0, 4(1, 1) = 1, or' e (cr,V^) = V^ ■ 

The last equality is related to the Dykhne theorem and 
generates a loop on the surface where (0) is satisfied. This 
loop is common for all such surfaces, i.e. they intersect 
each other on it. 

We computed 4(4 , 4) for the structures mentioned 
above: He, Fl, Co, Rh, and then extracted the parameter 
A corresponding to those structures. 

Before discussing these results we describe briefly the 
numerical algorithm which was used. This algorithm 
deals with hexagonal Bravais lattices represented as a 
grid of equilateral triangles. A subdivision procedure 
produces a set of N x N small similar triangles, the cen- 
ter of each one of them connected with the centers of 
3 nearest neighboring cells by simple resistors. The re- 
sistor connecting cells i and j has a resistance equal to 
(fjj + <jj)/(2*s/io-i<jj), where Ui is the conductivity of the 
cell i. Translational symmetry of the composite is re- 
flected by imposing periodic boundary conditions for the 
currents. The algorithm solves up to 10 6 linear equa- 
tions arising for the subdivided 3-color elementary cell. 
Our computational procedure always gives a sequence of 
effective conductivities er^, which converges monotoni- 
cally (as N — > 1000) from below. This fact was es- 
tablished by treating the solvable case following from 
Dykhne's theorem, where o~ e is known exactly. To get 
a sequence of upper bounds we then used the duality 
property, i.e. simulating the dual problem with conduc- 
tivities 1/4 1 1/4 an d consequently obtaining a mono- 
tonic convergence from above for the sequence <t^. 

The simulations were done on the 4 characteristic 
curves |^| on the surface where it intersects with the planes: 

1)4 = 0.001, 2)4 = 1, 3)4 = 4, 4)4 + 4 = 1. 

(16) 

These curves reflect the behavior of the surface relatively 
well in accordance with the chosen type of composite. 
From the assumption of the algebraic nature ([ll]) of those 
curves we were able to extract, for every type of compos- 
ite, a corresponding parameter A by the following pro- 
cedure. In a wide range of A, for each calculated point 
(4 7 i = 2,3) on the plane, relative deviations e were 
calculated 



3 Because of the divergence of the calculational procedure at 
a'i = 0, the first of the characteristic curves was calculated 
using (72 = 0.001, for which the results can still be trusted. 
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if cr+ < a e (A,a^) 



e(A,o<) = { if a N < <x e (A,ai) < a+ 

°\^ d if a e(A, CT-) < (Tjy 



(17) 



and a maximal value of those deviations £ m ax(A) was 



determined for every value of A by scanning over the 
entire area 0< a' 2 3 <1 . We then determined the best 
value of A by minimizing the function e max (A) . As shown 
in Fig. ||, A for the He and Fl structures is determined 
by very sharp minima. The other two minima are not as 
sharp. 



0. 1000 — 



o.ooio -J e 




FIG. 3. The fitting curves e ma x(A) for 4 different 
structures. The best value of A is determined by 
the minimum in each curve. Those values are: He - 
A He = 11.37, e£L=- 00026 ; Fl - A Fl = 
; Co - Aco = .305, eg° ax = m2 ; Rh - A Rh 



t Fl -- 

^max 



As one could expect, the Rh structure has A = since 
this is the unique value of A for which the solution of 
equation ([ll]) has the features typical of the structures 
with traps. The values of A which minimize e m ax(A) for 
the other structures are listed in the caption of Fig. |^. 

Figs. show the computed results for upper and 

lower bounds on the relative bulk effective conductivity 
a' e , i.e. <7jv(<T2, 03) and ajj (a' 2 , a' 3 ) at maximal N, for the 
4 microstructures of Fig. [l[ Note that these upper and 
lower bounds often appear to coincide due to insufficient 
resolution in the figures. The axes in the figures are la- 
beled with S.ef for the relative bulk effective conductivity 
<7g and si, s2, s3 for the component relative conductivities 
a{ = 1,02, a' 3 respectively. 

The pairs of almost merged points shown in Fig. 
correspond to barely separated upper and lower bounds. 




0.2 0.4 0.6 

Sl=l,s2=0,s3, (A=11.37) 



.2 0.4 0.6 
sW, s2=s3, (A=11.37) 



FIG. 4. Numerical results <r^(o"2, 03), 0jv(02i ct 3) (black 
dots) and analytic solution <j' e (o' 2 , o' 3 ) of cubic equation jll| ) 
with ^4=11.37 (bold line) at the four sections ( |l6| ) of the sur- 
face corresponding to the He - structure. 




0.4 0,6 0.8 1 
Sl=l,s2=0.001,s3, (A=3.76 



FIG. 5. Numerical results 0^(02,03)1 0^(02,03) (black 
dots) and analytic solution 0-^(0-2,03) °f cubic equation ( pl| ) 
with ^4=3.76 (bold line) at the four sections ( |l6| ) of the surface 
corresponding to the Fl - structure. 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 
sl'sM,s3, (A=.305) sW,s2=l-s3, (A=. 305) 




0.2 0.4 0.6 0.8 1 o 0.2 0.4 0.6 0.8 1 
sM,s2=0.001,s3, (A=.305) sW , s2 = s3 , (A=.305) 



FIG. 6. Numerical results o"jv( a 2i °"3)i °"jv(°2! 03) (black 
dots) and analytic solution o' e (a' 2 , o~' 3 ) of cubic equation 
with ^4=0. 305 (bold line) at the four sections of the sur- 
face corresponding to the Co - structure. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 
sl=s2=l, s3, (A=0) sM, s2=l-s3, (4=0) 




0.2 0,4 0,6 0,8 1 
sl=l, s2=s3, (A=0) 



FIG. 7. Numerical results o~^(o~2, 03), o~^,(a2,o-' a ) (black 
dots) and analytic solution a' e (a' 2 ,a' 3 ) of cubic equation ( p"l| ) 
with A=0 (bold line) at the three sections ( |l6| ) of the surface 
corresponding to the Rh - structure. The fourth dependence 
a' e (0,(73) =0 is omitted because of its triviality. 



V. EFFECTIVE CONDUCTIVITY OF RANDOM 
STRUCTURES: EXTENSION OF THE 
ALGEBRAICITY CONJECTURE 

In this chapter we discuss briefly a possible extension of 
the algebraicity conjecture for 2D three-component struc- 
tures, which are non-regular but macroscopically homo- 
geneous. On length scales large compared to the inhomo- 
geneities, we can characterize the macroscopic response 
of such a medium by a single number, the plane effective 
conductivity <r e . A reexamination of the basic properties 
(EH]) of a e (ai, 02, 03) shows that one of them (||) must 
be discarded, while the other three (homogeneity of 1-st 
order, duality, and compatibility) continue to be valid. 
It is important to note that these properties hold irre- 
spective of whether the microstructure of the (isotropic) 
composite is ordered or disordered. 

Instead of Eq. (Q), which implies full symmetry under 
the group S3 of all 3-color permutations, we first consider 
the case where the structure is only symmetric under 
the group C3 of cyclic 3-color permutations: 1^2^ 
3 — > 1. An example of such a microstructure is shown 
in Fig. ||. (Note the arrangement of "flowers" near the 
center of that figure.) It follows that <7 e ((7i, 02, 03) is a 
cyclic permutation invariant function. The group C3 is 
a subgroup of index 2 of the full permutation group S3, 
which characterizes all the regular structures He, Fl, Co, 
Rh discussed in the previous sections. 




FIG. 8. 3-color 2D composite of Fl structure locally dis- 
turbed in the centers of only three plaquettes with differ- 
ent colors. Such color tessellation preserves an equal fraction 
of every color and is invariant under the cyclic permutation 
group C3, but not under the full S3 permutation group. 

Being non-reflecting, the group C3 has more basic in- 
variants than the group S3 [see Eq. (fT(i|)l, namely 

h(C 3 ) = (Tl + (7 2 + CT3 , h(C 3 ) = (71(72 + O2O3 + (7 3 (7l , 
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h(Cs) = a x a 2 a-i , h(C 3 ) = {p\ - a 2 ){a 2 - cr 3 )(cr 3 - 01) . 

Nevertheless, the additional cubic invariant Ia{C^) can- 
not be incorporated into a cubic equation for <r e , because 
that would violate the duality requirement. Thus, we are 
lead back to the cubic equation ( pd|) , which is dictated 
not only by the strong requirement (W) of the full permu- 
tation invariance S3 , but even by the milder requirement 
of cyclic permutation invariance C3. 

In the case of a random composite, we usually char- 
acterize the microstructure by a statistical distribution 
function of the local conductivity, which can be either 
(j 1 , (72, or (73 at any point. Such a description results in 
an ensemble of representative structures, each one with 
its own form for the function a e (ai, o 2l (73). If we assume 
that the distribution function has the permutation sym- 
metry of either S3 or C 3 , then the ensemble average of 
<7 e ((7i, (72, cr 3 ) will have that symmetry too, even though 
individual samples may violate it. It follows that, the 
same considerations which led us to stipulate the form 
(|TT|), for the minimal polynomial equation that <7 e could 
satisfy for regular structures, also lead to that same equa- 
tion for random structures, if the statistical model for 
those structures is invariant under either S3 or C3. Nu- 
merical tests of this conjecture remain to be performed. 

VI. CONCLUSION 

In the present paper we have introduced the algebraic- 
ity conjecture for the effective conductivity problem of 
isotropic 2D three-component regular composites. This 
conjecture is based on the general properties which are 
satisfied by the effective conductivity a e (ai, o 2 , (73). The 
algebraic equation of minimal order for a e is a cubic equa- 
tion with 1 positive free parameter A responsible for the 
non-universality. This equation satisfies Dykhne's the- 
orem (0) independently of A and has only one positive 
root, thus avoiding the possibility of multiple physical 
solutions. The value A — 1/3 corresponds exactly to 
the Bruggeman effective medium approximation for a 2D 
composite with 3 equally partitioned components. 

We have found support for this conjecture by numeri- 
cal calculations on four different infinite 2D 3- color class 
equivalent regular structures of P6mm(L)|P6mm(L ) and 
P6(L)|P6(L') types: He, Fl, Co, Rh. We have established 
that cr e (<7i, (72, (73) is a non- universal function with essen- 
tial dependence on the microstructure even for totally 
symmetric structures: 

1) The cubic equation ( pi] ) with A=11.37 governs the 
conductivity problem in He structure with a very high 
precision e max ~ 1CT 4 . 

2) There is good agreement (e max ~ 10 -3 ) between the 
cubic equation with A—3.1Q and numerical results for the 
Fl structure. 



3) In the Co structure the estimated value A=.305 is 
near 1/3, which would follow from the Bruggeman ef- 
fective medium approximation. This may indicate some 
similarity between the conducting properties of a 3- 
component random microstructure and those of the or- 
dered Co structure. 

4) The Rh structure needs special attention. It belongs 
to the family of structures with unicolor traps, i.e. with 
structures where the presence of just one non-conducting 
((7i = 0) component is enough to make the composite 
an insulator (o~ e = 0). In percolation theory this corre- 
sponds to a threshold P3 = 1/3, in contrast with a 2-color 
composite where p 2 — 1/2. If the cubic equation is valid 
for this structure with A=0, then a e = a\a 2 a3. Unfor- 
tunately, these computations are unable to resolve this 
question for the Rh microstructure. 

5) The different microstructures which belong 
to a common plane symmetry group (He, Rh - 
P6mm(L)|P6mm(L'), Fl, Co - P6(L)|P6(L')) are charac- 
terized by distinct values of A. This means that this pa- 
rameter has a topological nature and is sensitive to more 
than just the symmetry properties of the elementary cell 
or the plane group of the entire color lattice. 

Finally we discussed a possible extension of the alge- 
braicity conjecture to other types of 2D three-component 
structures. First we showed that even if the microstruc- 
ture is only invariant under the C3 subgroup of the S3 
permutation group, a cubic equation for a e must still 
have the form ([Ti|). We then extended this conjecture 
also to the case of random microstructures, described by 
any statistical model that is invariant under either S3 or 
C3. To test this idea, it would be useful to do numerical 
calculations on such structures. This is left for future 
investigations. 
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APPENDIX A: 

We now give a simple proof of the duality relations for 
a 2D isotropic medium composed of an arbitrary number 
n of isotropic components. 

Suppose a 2D medium with a continuous distribution 
of conductivity <r(r) is subjected to an average electric 
field E e . The system of equations consists of Ohm' law 



J(r) = a(r) E(r) , 
and the equations for local fields 

V x E(r) = , V ■ J(r) = 



(Al) 



(A2) 



with appropriate boundary conditions on the electrical 
potential. 

We are interested in the relation between the current 
J e averaged over the system, J e = S^ 1 J J(r)dS and 
the averaged field E e = J~E(r)dS. By virtue of the 
linearity of (Al, A2) this relation will also be linear 



(j P E P 



CT e {CT(r)} 



(A3) 



where the tensor of the bulk effective conductivity a e — 
a e {cr(r)} is actually a tensorial functional. In the case 
of n homogeneous anisotropic components, the latter be- 
comes a tensorial function a e (a\, a 2 , <7„). Further 



simplification arises when all components are isotropic — 
<7 e ((71,02, ...,er„) and, finally when the entire composite 
is also an isotropic medium — a e (a±, cr 2 , ...,a n ). 

In order to transform to the dual problem we rotate 
the x, y -components of J and E by 90° in the plane 

J'(r) = .RE(r), E'(r) = )?J(r) , J'(r) = a'(r) E'(r) , 



ff'(r) =Ra~ 1 (r) R^ 1 , R = 



-1 

1 



(A4) 



Eqs. A2 are thereby transformed as follows 

V • J'(r) = V x E(r) = , V x E'(r) = V • J(r) = , 

while Jg and Eg are connected by 

J e = CT e Eg , d' e =R d^R- 1 , (A5) 
where the bulk effective conductivity tensor of the dual 



problem a' e is defined in accordance with (A3) 

d' e = d e {d'(v)} = d e {r d-\r)R- 1 } . (A6) 

The last two equations ( A5, A6|) give a duality relation 

CT e {Ra^i^RT 1 } -R-a e {(?(r)} • R 1 = I, (A7) 

where / is the unit matrix. If, instead of the continuous 
fields tr(r) we deal with anisotropic media composed of n 
homogeneous anisotropic components, then 

(jg (R fff 1 i? (j^ 1 R- 1 ,...,^,^- 1 R- 1 ) ■ R ■ 

■&e (04,(72, . ..,(?„) ■ R^ 1 = I . (A8) 

In the case of an anisotropic composite with n homoge- 
neous isotropic components we will have 

<7 e (erf 1 , <7,T , er^ 1 ) • R ■ a e {<j\, a 2 , cr„) • R^ 1 = I . 

(A9) 

The principal values of a e satisfy Keller's theorem 
5ffa:W,.",0 ■dt v (a 1 ,a 2 ,...,a n ) = 1 , (A10) 

In general, the directions of the principal axes depend 
on the values of a\,a 2 , etc. But when the symmetry of 
the microstructure is sufficiently high, those directions 
will be fixed by that symmetry. Finally, in the case of 
an isotropic composite, the last equations reduce to the 
self-duality relation 

Ozip^ 1 ,a 2 X , .:,cr^ 1 ) ■ a e {ai,a 2 , ...,cr n ) = 1 . (All) 
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